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Direct observation of magnetic cycles in Sun-like stars 
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A sample of 19 solar-type stars, probing masses between 0.6 and 1.4 solar mass and rotation periods between 3.4 and 
43 days, was regularly observed using the NARVAL spectropolarimeter at Telescope Bernard Lyot (Pic du Midi, France) 
between 2007 and 2011. The Zeeman-Doppler Imaging technique is employed to reconstruct the large-scale photospheric 
magnetic field structure of the targets and investigate its long-term temporal evolution. We present here the first results of 
this project with the observation of short magnetic cycles in several stars, showing up a succession of polarity reversals 
over the timespan of our monitoring. Preliminary trends suggest that short cycles are more frequent for stellar periods 
below a dozen days and for stellar masses above about one solar mass. The cycles lengths unveiled by the direct tracking 
of polarity switches are significantly shorter than those derived from previous studies based on chromospheric activity 
monitoring, suggesting the coexistence of several magnetic timescales in a same star. 
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1 Introduction 

Sun-like stars are characterized by convective envelopes in 
which large-scale plasma flows (related, in particular, to ra- 
dial and latitudinal differential rotation and to the Coriolis 
force) are able to trigger a global dynamo (Parker 1955). 
This continuous generation of a large-scale field is related 
to surface variability affecting a wide range of temporal and 
spatial scales, including quasi-periodic polarity reversals as- 
sociated to magnetic cycles. Recent numerical models, in 
particular global MHD simulations, are able to mimic some 
characteristics of this cyclic behaviour for Sun-like stars 
(Ghizaru et al. 2010, Brown et al. 201 1). 

This magnetic activity is observable through many ob- 
servational proxies, from the photosphere (e.g. broad-band 
visible photometry) to the corona (X-ray and radio emis- 
sions). The long-term monitoring of chromospheric activ- 
ity of tens of solar- type stars, carried out at Mount Wilson 
since 1965 (Wilson 1978) has allowed for the detection of 
periodic variations in a number of objects (Baliunas et al. 
1995). Cyclic patterns were observed with a variety of cy- 
cle lengths, as well as seemingly more erratic activity fluc- 
tuations in other objects (in particular young stars) or no 
detectable activity at all in some others. More recently, as- 
teroseismology was demonstrated its ability to investigate 
magnetic variability through variations of the p-mode am- 
plitudes and frequencies (Garcia et al. 2010) with the detec- 
tion of variations of these quantities over a few tens of days, 
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presumably linked to magnetic activity, for the F5 V star HD 
49933. 



In addition to this wealth of indirect indicators of stel- 
lar activity, spectropolarimetry now enables us to perform 
direct measurements of surface magnetic fields and follow 
the long-term temporal evolution of large-scale magnetic 
geometries. So far, it allowed the observation in Sun-like 
stars of one global polarity switch (Petit et al. 2009) and of 
a full magnetic cycle (Fares et al. 2009). 



Our aim is to study the long-term variations of the mag- 
netic field properties of a sample of solar- type stars, using 
both direct and indirect measurements. Our observed sam- 
ple includes 19 FGK-type stars on the main sequence, mon- 
itored since 2007. We probe here stellar masses between 0.6 
and 1.4M , and rotation periods between 3.4 and 43 days. 



After a brief description of the instrumental setup, data 
reduction and multi-line extraction of Zeeman signatures, 
we explain the reconstruction technique of the large-scale 
topology of the stars, and the computing of a chromospheric 
activity indicator. We then highlight three representative ex- 
amples of different types of variability observed in our sam- 
ple. We finally discuss the results derived from our measure- 
ments. 
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2 Instrumental setup, data reduction, and 
extraction of Zeeman signatures 

We use data from the NARVAL spectropolarimeter (Auriere 
2003), installed at Telescope Bernard Lyofl (Pic du Midi, 
France). The instrumental setup is strictly identical to the 
one described by Petit et al. (2008). The spectrograph unit of 
NARVAL benefits from a spectral resolution of 65,000 and 
covers the whole wavelength domain from near-ultraviolet 
(370 nm) to near-infrared (1,000 nm). Thanks to the polari- 
metric module, NARVAL can provide intensity, circularly 
or linearly polarized spectra. In the present study, we restrict 
the measurements to Stokes I and V. 

The circularly polarized spectra allow the detection of 
large-scale photospheric magnetic fields, thanks to the Zee- 
man effect. However, when observing cool dwarfs, the sig- 
nal-to-noise ratio of circularly polarized spectra produced 
by NARVAL is not sufficiently high to reach the detection 
threshold of typical Zeeman signatures (which amplitude 
does not exceed 10 _4 / c for low-activity stars, where I c is 
the continuum intensity). To solve this problem, we calcu- 
late from the reduced spectrum a single, cross-correlated 
photospheric line profile using the Least-Squares-Decon- 
volution (LSD) multi-line technique (detailed by Donati et 
al. 1997 and Kochukhov et al. 2010). Thanks to the large 
number of available photospheric lines in cool stars (several 
thousands in the spectral domain of NARVAL), the noise 
level is reduced by a factor of about 30 with respect to the 
initial spectrum. As an illustration, Fig. \T\ shows the result- 
ing LSD signatures for successive observations of the Sun- 
like star £ Boo A. 

3 Magnetic mapping and chromospheric 
emission 

The Stokes I and Stokes V LSD profiles allow the derivation 
of various quantities to study the temporal variations of the 
magnetic field properties. Here we focus on the reconstruc- 
tion of the surface distribution of the magnetic vector and 
on the computation of a chromospheric activity index. 

3.1 Magnetic maps 

To reconstruct the surface magnetic geometry of the stars, 
we use Zeeman-Doppler Imaging (ZDI). This tomographic 
inversion technique is based on the modelling of the rota- 
tional modulation of the circularly polarized signal (Semel 
1989). The time series of polarized signatures are iteratively 
compared to artificial profiles corresponding to a synthetic 
magnetic geometry, until a good fit is obtained between the 
model and the observations (Donati & Brown 1997, Donati 
et al. 2006). Thus, ZDI enables to recover, to some extent, 
the location of magnetic regions, as well as the strength and 

1 The Bernard Lyot Telescope is operated by the Institut National des 
Sciences de l'Univers of the Centre National de la Recherche Scientifique 
of France. 
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Fig. 1 Normalized Stokes V profiles of £ Boo A for the 
summer of 2007, after correction of the mean radial velocity 
of the star. Black line represent the data and red lines corre- 
spond to synthetic profiles of our magnetic model. Succes- 
sive profiles are shifted vertically for display clarity. Rota- 
tional phases of observations are indicated in the right part 
of the plot and errors bars are illustrated on the left of each 
profile. 



orientation of the magnetic vector in magnetic spots. The 
application of this technique to cool stars with low v sin i 
and moderate to low magnetic activity is described by Petit 
et al. (2008). In this case, ZDI is only sensitive to low-order 
field components, contrary to the chromospheric flux which 
includes also the contribution of smaller scale magnetic el- 
ements. 

The resulting maps for the three stars presented here are 
illustrated in Fig.[2l[|andH 
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Table 2 Chromospheric activity indices for HD 78366, 
HD 190771 and £ Boo A for each corresponding set of ob- 
servations. 



Star 


Set of obs. 


NcallH 


HD 78366 


2008.09 
2010.04 
2011.08 


0.273 ± 0.004 
0.291 ± 0.006 
0.278 ± 0.003 


HD 190771 


2007.59 
2008.67 
2009.47 
2010.50 


0.335 ± 0.006 
0.338 ±0.011 
0.337 ± 0.007 
0.345 ± 0.006 


£ Boo A 


2007.59 
2008.09 
2010.48 
2010.59 


0.443 ± 0.008 
0.420 ± 0.008 
0.403 ± 0.006 
0.402 ±0.011 



3.2 NcaiiH -index 

From the Stokes I profiles, we construct an index to quan- 
tify the chromospheric emission changes in the Call H line. 
The complete pipeline of the computation is described in 
details in a forthcoming paper (Morgenthaler et al., submit- 
ted). We follow the methods of Duncan et al. (1991) and 
Wright et al. (2004), who define indexes based on Mount 
Wilson observations, and we calculate a Nc a iiH -index for 
our NARVAL observations. To help comparing our chromo- 
spheric estimate with older studies, we calibrated the index 
against the values derived at Mount Wilson. The Ncaiw- 
indexes we obtained for HD 78366, HD 190771 and f Boo 
A are detailed in Tab. [2 

4 Results 

Since the monitoring began a few years ago, long-term chan-| 
ges in the magnetic properties become observable in some 
of our targets. Both the magnetic quantities derived from 
ZDI and the chromospheric index exhibit temporal fluctua- 
tions over a wide range of timescales, due to rotational mod- 
ulation and longer-term magnetic trends. Three representa- 
tive examples of the different kinds of stellar variability we 
observed is described hereafter. 

4.1 Short magnetic cycle : HD 78366 

A simple type of variability is observed for HD 78366. This 
target is slightly more massive than the Sun and rotates about] 
two times faster (Tab.Q]). The data sets of this object are col- 
lected over three distant seasons. The corresponding mag- 
netic maps are shown in Fig. [2J We observe two polarity 
switches, especially visible in the polar area of the radial 
field component, which is of negative polarity in 2008.09, 
positive in 2010.04 (and associated at that time with a more 
complex magnetic field geometry), and negative again in 
201 1 .08. After the two observed polarity reversals, the mag- 
netic field retrieves its initial configuration. Assuming that 
the magnetic variability of this star is not much faster than 



the temporal sampling imposed by the right ascension of the 
star (visible only during winters), this first time- series sug- 
gests that HD 78366 may obey to a magnetic cycle of about 
three years. We note that the chromospheric activity indica- 
tor NcaiiH (Tab. [3 seems to increase with the complexity 
of the large-scale magnetic field, with a clear maximum in 
2010.04. 

4.2 Fast polarity reversals : HD 190771 

A more complex type of variability is illustrated by HD 
190771. It has a mass similar to the Sun's, but has a rota- 
tion period of 8.8 days (Tab.Q]). In Fig. [3] we plot the mag- 
netic maps derived for this star. A polarity reversal is visible 
on the strong azimuthal component between 2007.59 and 
2008.67 (Petit et al. 2009). Between 2008.67 and 2009.47, 
the magnetic geometry changed in a different manner : the 
magnetic field which was mainly toroidal in 2008.67 be- 
came mostly poloidal in 2009.47. A second polarity reversal 
took place between 2009.47 and 2010.50, this time on the 
radial field component. In this case, the two successive po- 
larity switches do not imply that the initial magnetic state is 
reached again, so that the observed variability is not taking 
the form of a cycle. 

In addition, we observe that the magnetic field intensity 
is correlated with the chromospheric emission. In the first 
three years, both the field strength and the chromospheric 
flux are roughly stable. We then observe a significant in- 
crease of these quantities in 2010.50 (Tab. [2]). 

4.3 Fast and complex variability : £ Bootis A 

Finally, another, more complex type of variability is ob- 
served with £ Boo A, the less massive and most rapidly ro- 
tating star of our three examples (Tab. [I]). It was observed 
at seven epochs, for which the magnetic field geometry and 
NcaiiH were derived (Morgenthaler et al., submitted). Here| 
we highlight two results of this long-term monitoring. 

The first one refers to the 2007.59 and 2008.09 data sets 
(top part of Fig. HI). We observe that within a six months in- 
terval, the intensity of the magnetic field decreased by about 
50% and that the magnetic geometry, which was quite sim- 
ple in 2007.59 with an aligned dipole and a prominent ring 
of azimuthal field, became more complex and less axisym- 
metric in 2008.09, with a less pronounced toroidal surface 
component. The decrease of the large-scale field strength is 
also observed in the chromospheric flux, with a sharp drop 
of emission between the two epochs (Tab. [3. 

The second example is visible in the set of observations 
collected during the summer of 2010, which we decided to 
split in two subsets (2010.48 and 2010.59) to take into ac- 
count the fast variations of the Zeeman signatures over this 
short timespan. In the correponding magnetic maps (bottom 
part of Fig. [4]), the most striking evolution is a sharp in- 
crease of the azimuthal magnetic field. These changes are 
taking place at a roughly constant level of chromospheric 
emission. 
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rot. phase rot. phase rot. phase 



Fig. 2 Magnetic maps of HD 78366, derived from 2008.09, 2010.04 and 201 1.08 observations (from left to right). For 
each data set, the 3 charts illustrate the field projection onto one axis of the spherical coordinate frame with, from top to 
bottom, the radial, azimuthal, and meridional field components. The magnetic field strength is expressed in Gauss. 
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Table 1 Fundamental parameters of HD 78366, HD 190771 and £ Boo A. T e ff, mass, radius and v sin i are taken from 
Valenti & Fischer (2005), Fernandes et al. (1998), Gray (1984), Petit et al. (2005). The rotation periods and inclination 
angles are derived from ZDI and from Toner & Gray (1988). 



Star 


T eff 


Mass 


Radius 


vsirn 


peg 

r rot 


inclination 




(K) 


(M ) 


(Rq) 


(km/s) 


(d) 


(degrees) 


HD 78366 


6014 ± 50 


1.34 ±0.13 


1.03 ±0.02 


3.9 ±0.5 


11.4 ±0.1 


60 ± 15 


HD 190771 


5834 ± 50 


0.96 ±0.13 


0.98 ± 0.02 


4.3 ±0.5 


8.8 ±0.1 


50 ± 10 


£ Boo A 


5551 ± 20 


0.86 ± 0.07 


0.80 ± 0.03 


3.0 ±0.4 


6.43 


28 ±5 




£ Boo A is therefore submitted to fast and complex sur- 
face changes that are different from those of the two previ- 
ous stars, and reminiscent of the complex behaviour of other 
rapid rotators observed in the past (e.g. Donati et al. 2003). 

5 Discussion 

All stars of our sample show variability over the four years 
of our monitoring, but of different types. Stars which show 
at least one field reversal over this timespan have in com- 
mon a fast rotation period (at least twice the solar one) and 
masses equal or slightly larger than that of the Sun (Fig. [5]). 
In Fig. [51 we include r Bootis, which is not part of our sam- 
ple but which is reported to be affected by a short magnetic 
cycle of two years at most (Fares et al. 2009). We stress also 
that active stars with masses below our lower mass bound- 
ary (in particular, mid-M dwarfs with masses just below the 



fully convective limit) are reported to possess strong, simple 
and stable surface magnetic fields (Morin et al. 2008a,b). 

r Boo and HD 78366 were also observed at Mount Wil- 
son as chromospherically active stars. For r Boo, Baliunas 
et al. (1995) report a cycle of twelve years, versus two years 
from spectropolarimetry. Concerning HD 78366, periods of 
six and twelve years were identified using the Mount Wil- 
son time- series, against about three years in our investiga- 
tion. We therefore note that, at least for these two examples, 
the cycle lengths derived by chromospheric activity seem 
to be longer than those derived by spectropolarimetry. We 
suggest that this apparent discrepancy may be linked to the 
different temporal sampling inherent to the two approaches, 
so that the sampling adopted at Mount Wilson may not be 
sufficiently tight to unveil short activity cycles. 

Future observations of our stellar sample will allow us 
to investigate longer timescales of the stellar magnetic evo- 
lution. The sample includes several solar twins (Petit et al. 
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Fig. 5 Rotation period versus mass for the stellar sam- 
ple. Pink symbols stand for stars with at least one polarity 
switch. 



2008) which have not showed cycles yet, and which will 
help us to determine how small departures from the solar 
fundamental parameters may affect the characteristics of 
magnetic cycles. More generally, a regular monitoring of 
our targets over more than one decade will enable us to de- 
termine more precisely the relation between the length/oc- 
currence of magnetic cycles and the rotation/mass of Sun- 
like stars. 
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